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Abstract 

We dynamically realize the super symmetric inverse seesaw mechanism in the next to the minimal 
supersymmetic standard model (NMSSM) in a minimal form. The crucial observation is that the 
inclusion of a dimension-five operator N 2 S 2 / 'M* , which violates lepton number and is suppressed 
by the fundamental scale M*, can readily produce a light neutrino mass scale m u ~ 0.01 eV. 
This model is protected by Z± x Z^ 1 discrete symmetry and a very predictive parameter space is 
restricted by the requirement of naturalness. Then we study some interesting phenomenological 
consequencese: (i) A light sneutrino (~ 8 GeV), which is the lightest supersymmetric particle 
(LSP) and then the dark matter candidate, can elegantly explain the DAMA/CoGeNT results and 
even account for the FERMI-LAT gamma ray signals; (ii) It is possible that the sneutrino is an 
asymmetric dark matter and thus the SIdm : ~ 5 : 1 naturally predicts a light dark matter; 
(iii) Higgs physics maybe dramatically changed due to the significant coupling y^LH u N and the 
light sneutrinos out of the inverse seesaw sector, consequently the lightest Higgs boson is likely 
dominantly decay to a pair of light sneutrinos. In the DAMA/CoGeNT region, the Higgs typically 
decays to a pair of LSPs with branching ratio > 90%. 

PACS numbers: 12.60.Jv, 14.70.Pw, 95.35.+d 
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I. INTRODUCTION AND MOTIVATION 



Theoretically, imposing supersymmetry on the standard model (SM) leads to the mini- 
mal supersymmeric standard model (MSSM), which is the most promising way to stabilize 
the weak scale Mz ~ 100 GeV <C Mpi. But phenomenologically the MSSM suffers the 
[i— problem and the little hierarchy problem. As a popular solution of these problems, the 
next to the minimal supersymmetric standard model (NMSSM) [1] has recently attracted 
much attention. Its superpotential and soft terms for the Higgs sector are given by 

Wnmssm =^SH u ■ H d + —S , (1) 
-C soft =m 2 Hu \H u \ 2 + m 2 Hd \H d \ 2 + m 2 s \S\ 2 + (\A X SH U ■ H d + ^A K S 3 + h.c^j , (2) 

where H u ■ H d = H\H\ — H 2 H\ (for simplicity the dot will be omitted hereafter). Once 
the singlet Higgs field develops a non-zero vacuum expectation value (VEV) v s = (S) and 
breaks the Z 3 symmetry around the TeV scale, the fi term is generated as fi — Xv s , which 
gives an elegantly dynamical solution to the /i— problem (and also the fi/B/j, problem). 

However, the NMSSM cannot possibly be the whole story because a tiny non-zero neutrino 
mass < 0(0.1) eV is needed to be generated. Classically, such a tiny mass parameter can 
be addressed by the dimension-five operators (H u Li)(H u Lj) / A, with A ~ 10 14 GeV. To 
understand the UV origin of these high-dimensional operators, the elegant seesaw mechanism 
is proposed [2]: 

M ij 

^— = isw, + -^w. (3) 

In this mechanism, three right-handed neutrinos (RHNs) with large Majorana masses ilf^ ~ 
A (the seesaw scale) are introduced. Although this mechanism is attractive for its simplicity 
and well motivated in the 5*0(10) grand unification theory (GUT), generically speaking it is 
not testable at the future colliders. If the seesaw scale is lowered to TeV, the RHNs will have 
very weak Yukawa couplings y 1 ^ ~ 10 -6 and thus it is still irrelevant for collider experiments 
(note, however, in the NMSSM with SN 2 coupling to account for the TeV-scale seesaw, the 
RHNs may be accessible [3]). Moreover, the origin of such an intermediate seesaw scale 
remains a puzzle. 

The inverse seesaw [4] can lower the seesaw scale while keep the large Yukawa couplings 
for the left-handed neutrinos. The superpotential is given by 

W IVS = y N H u LN + M N NN + ^N 2 . (4) 
Here, for simplicity, we neglect the family indices. The light neutrinos are the three lightest 
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eigensates of the 9x9 mass matrix, which, in the basis (ul, N, N), is given by 

/ m T D \ 



M 



inverse 



m D M% 
\ M N M R J 



(5) 



where each block denotes a 3 x 3 matrix with Dirac neutrino mass matrix mp = Y^(H^). 
In this paper, to illustrate our main idea, we consider one family of neutrinos for simplicity. 
Then the effective light neutrino, which is dominated by u^, obtains a mass 

We are interested in the case with at the weak scale and ~ 0(1). So, the smallness of 
neutrino masses are generated by small Majorana masses M^. Since this mass term violates 
lepton number, we see that the small neutrino masses may be ascribed to the degree of 
lepton number violation. 

However, such an inverse seesaw mechanism suffers from a fine-tuning problem, namely 
the Majorana mass term must be very small. Some efforts have been tried to construct 
dynamical models to solve this naturalness problem. Commonly, they start from a lepton 
number (L) conserving model and then the smallness of Mfj is a result of the spontaneous 
breaking of U{1)l [5]. Usually, some new complicated structures and degrees of freedom are 
introduced to realize the dynamical goal. In this work we start from a L— violating theory 
and present a very simple model to realize the inverse seesaw without introducing any new 
extra fields or any new ad hoc terms. Thus, our framework is a kind of minimal dynamical 
realization for the inverse seesaw mechanism. 

Our model has rich and interesting phenomenology. First of all, a left-right mixing 
sneutrino as the lightest supersymmetric particle (LSP) may serve as a dark matter (DM) 
candidate. In fact, the recent DM direct detection results from DAMA [6] and CoGeNT 
[7] can be explained by a light DM with a mass around ~ 8 GeV and a rather large DM- 
nucleon scattering cross section asi ~ 10 _40 cm 2 [8] (a recent interesting work on light DM 
with isospin violation can recoil DAMA/CoGeNT with other experiments is explored in the 
[9]). Such a light DM candidate is absent in the MSSM [10], while in the NMSSM its presence 
requires some fine-tuning in the parameter space [11] (some other exotic extensions allowing 
for a light DM can be found in [12]). However, in our model, if the light sneutrino LSP has 
a proper doublet sneutrino component, it can possibly explain the DAMA/CoGeNT signals 
by the Z— boson mediated recoil [13]. In some cases, it may have leptonic annihilation into 
t + t~ with a rate ~ 10~ 27 cm 3 /s, which, according to the analysis in [14], may be favored 
by the FGST (FERMI Gamma-Ray Space Telescope) observation of gamma ray in the 
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region 0.3 — 100 GeV (please see some following study in [15]). On the other hand, since 
the sneutrino LSP carries lepton number, it can be the asymmetric dark matter (ADM) 
candidate [16-18]. Finally, provided that the sneutrino is so light, the Higgs boson may 
mainly decay invisibly into a pair of light sneutrinos. Such an interesting phenomenology 
not only affects the collider physics about Higgs but also may provide a way to alleviate the 
little hierarchy problem. And in the DAMA/CoGeNT region, the Higgs typically decays to 
a pair of LSPs with branching ratio > 90%. All these phenomenological consequences will 
be discussed in our analysis. 

This work is organized as follows. In Section II a dynamical model for the inverse seesaw 
is constructed and some discussions about the naturalness of the model are presented. In 
Section III, we explore some interesting phenomenology for this model. The conclusion 
is given in Section IV. Some necessary interactions of the sneutrinos are presented in the 
appendix. 

II. MINIMAL DYNAMICAL MODEL FOR THE INVERSE SEESAW 
A. The Model and Its Symmetry 

As stated in the preceding section, the inverse seesaw mechanism suffers from the fine- 
tuning problem. To avoid this problem, we now construct a simple dynamical model based 
on the NMSSM. First, we see from Eq. (6) that with the assumption 3> M^, the light 
neutrino mass is then an absolute result of the smallness of Mfj ~ (9(0.05) eV. In particle 
physics, it is very illustrative to get such a small mass scale from some 'fundamental' mass 
scale, namely 



where M* is the fundamental scale such as the Planck scale or the string scale. In the usual 
GUTs, the effective fundamental scale can be decrease by a factor about 3 due to the large 
degrees of freedom at the GUT scale [19]. So M* can be around 10 17 GeV in principle. 
Therefore, it naturally guides us to conjecture that M# may simply open a window to the 
fundamental scale (UV) physics. In this spirit, we construct the following simple model to 
dynamically realize the inverse seesaw mechanism 



(5 TeV) 2 



0.01 eV , 



(7) 



-C 



w DIS 




(y N LH u N + X 1 SNN) + -^fS 2 N\ 

m 2 L \L\ 2 + m 2 N \N\ 2 + m%\ft\ 2 + (y N A N H u LN + X^SNR + h.c. 



) 



(8) 



(9) 
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Such a very simple model can provide all the elements we need and it must be the minimal 
framework to dynamically realize the inverse seesaw mechanism because we do not introduce 
any new fields (only the NMSSM plus right-handed neutrinos). With the singlet S developing 
a TeV-scale VEV v s , both M N and Mfj can be generated 

M N = X lVs ~ 100 GeV , (10) 

Mr = l^T ~ °- 01 eV ' 

where the operator coefficient A2 is set to be 1. So, can be naturally at weak scale and 
Mfi can be close to the light neutrino mass scale. 

Now, the fermion spectrum contains a light Majorana neutrino with mass given by Eq. (6), 
whose corresponding eigenstate is 

v\ ps sm9iv L — cos#iiV, (12) 



with sin Q\ ps Mn / ym^+M^. The other two orthogonal Weyl fermions form a Dirac 
fermion: 

1/2,3 ~ —/= (cosO^l ±iV + sin#iV) , (13) 
v2 



M 2 = yjM% + m 2 D + 0{Mff) PS M N . (14) 

We make a comment on the dimension-five operator introduced in our model. In fact, this 
operator closely resembles the classical dimension-five operator in the SM or MSSM: 

(LH U )(LH U ) NNSS 

— a y ^r- (15) 

The canonical seesaw scale is a somewhat peculiar scale A ~ 10 14 GeV because of the 
VEV v u < 174 GeV. But by virtue of the inverse seesaw mechanism, the smallness of 
neutrino mass is transmitted to the small Majorana mass of the extra singlet N. Combined 
with the NMSSM singlet Higgs field S that has a TeV-scale VEV, then we can obtain 
the neutrino mass scale naturally via the dimension-five operator introduced in our model 
without invoking new unknown scales. 

We make a comment here. The assumption mo > Mn is not realistic since it will lead 
to the light neutrino dominated by the TV component while the major doublet component 
as heavy as Mjv- As a result, electroweak phenomenology such as W, Z boson decay will 
be inconsistent with the current experimental data. On the other hand, if M^jm^, is too 
large, M$ will not be sufficiently generated in our scenario. So in the actual calculation, 
Mfq ~ 5m o will be assumed. 



5 



Note that our dynamical scheme is highly constrained. First, we need the NMSSM with 
v s ^> v where v = a/ (H®) 2 + (-^d) 2 - ^ * s dear that the larger M^v, the smaller m u . In 
the fully inverse seesaw limit, v s takes a value around 5 TeV to produce a neutrino mass 
scale m v ~ 0.05 eV which is the heaviest neutrino mass observed by the neutrino oscillation 
[20]. As a result, if the little hierarchy problem does not appear, we have fj, = Xv s ~ 100 
GeV with A ~ 10~ 2 . Further, we are interested in Mjv is several hundred GeVs, so typically 
Xi < 10 _1 and in most cases we shall ignore the role of this parameter. 

Since the small neutrino mass scale is simply a relic of fundamental scale physics, this 
inspires us to investigate the possible model at the fundamental scale. We find that a 
new U(1)' R symmetry can guarantee the general form of our model. At the renormalizable 
level, the supersymmetric model described by Eq. (1) plus Eq. (8) possess an accidental 
U(1)b x U(1)l x U(1)r symmetry with the field charges assigned as 

L : H u [0], H d [0], S[0], L[l], E°[-l], N[-l], N[l], 

B : H u [0], H d [0], S[0], Q[l], U c [-1], D c [-1], N[0], 

R: H u [2/3], H d [2/3], S[2/3], L[2/3], E c [2/3},--- (16) 

where the dots denote that all other fields carry the same charge 2/3. In fact, the U(1) R 
charge assignment is not fixed according to this superpotential and in the above we simply 
choose one as an example, which is consistent with SU(5) GUT. Note that the Z 3 symmetry 
simply is an accidental result of U(1) R symmetry, which forbids the bare mass terms. At 
the dimension- five level, the operator S 2 N 2 violates the global symmetry U(1)l and U(1) R 
simultaneously, but still leaves a discrete C U(1)l and a new U(1)' R invariance: 

R'^R- l -L. (17) 

If U(l )' R is generated to all orders, then it is not difficult to prove that one of its consequence 
is that U(1)b and the matter parity Z^ 1 = (— 1) 3 ( B_L ) conserve to all orders [21]. 

Some comments are in order. First, we only utilize the universal Abelian U(1)l charge 
invariance at the renormalizable model. But R— charge does not re-scale, since we fix it to be 
2 for the superpotential, not zero, which is quite different from the ordinary U(l). Second, 
we may directly use the U(1) R alone to accomplish this task, namely to pick out some 
required operators and forbid the others, rather than using such a redefined U(1) R > ( This 
means the basic symmetry of the whole model is only U(1) R x G gauge . This is an interesting 
conjecture, even if taking supersymmetry (SUSY) breaking into account, where the U{1) R is 
required [22]. Of course, according to the well known conclusion that no continuous global 
symmetry is a fundament symmetry, we may finally break it by higher dimension operators) . 
Finally, our model does not conserve the U{1)l and thus the U(1)b-l, which is an essential 
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difference from the dynamical inverse seesaw model proposed in the literature [5] . Although 
the U(1)b-l is n °t conserved, the matter parity is still well defined, which is from the fact 
Z2 (equally Z;f _i ) is left. We shall see below that this symmetry is important when it is 
regarded as an exact symmetry even if taking non-renormalizable operators into account. 

B. Tadpole Term for S 

A light gauge singlet at low energy may destabilize the weak scale. Without any fun- 
damental symmetry like gauge symmetry to protect the singlet, some dangerous tadpole 
terms with quadratic divergence may be generated by the Planck-scale-suppressed opera- 
tors [23]. Although in this section we consider the Planck scale suppressed operators as in 
the literatures, the results hold if we change the Planck scale to the effective fundamental 
scale. 

We first briefly review how the NMSSM can avoid the dangerous tadpole and mean- 
while generate a safe tadpole of S to solve the domain-wall problem associated with the 
spontaneous breaking of the low-energy accidental Z 3 discrete symmetry. In particular, 
the odd-dimensional operators in the superpotential and even-dimensional operators in the 
Kahler potential will not destabilize the electroweak scale [24]. 

The study in [25] utilizes a Z^ 1 x Zf symmetry as the fundamental symmetry of the 
NMSSM superpotential and Kahler potential even at the Planck scale. The Z± symmetry 
is a discrete subgroup of U(1)r, under which all NMSSM superfields transform as $ — > i$ 
and the superpotential transforms as W — > — iW. This symmetry forbids any terms with 
odd powers of superfields in the Kahler potential and even powers of superfields in the 
superpotential. When coupling the theory to the supergravity, the Z^ symmetry will be 
broken by SUSY-breaking effect and consequently singlet S will acquire a tadpole fi 2 S 
induced by the non-renormalizable operators. However, according to the conclusion in [24], 
such a tadpole is safe and able to solve the domain-wall problem since n 2 s S breaks Z 3 
explicitly at a proper level. 

In our work, the situation is different in two aspects. On the one hand, the new singlets 
N and N are charged under the matter parity Zif , which is regarded as an exact symmetry 
in our model (no R— parity breaking operators are introduced) and thus the corresponding 
tadpoles are forbidden. In other words, Z^ 1 protects N and N from tadpole. On the other 
hand, one essential point of our model is the operator S 2 N 2 /Mp\ with even powers of super- 
fields in the superpotential. So, we have to check whether it will induce a harmful tadpole 
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of S. In fact, the Z± C U(1)r symmetry is still conserved, with the charge assignment as 



Then W — > — W under this transformation. In the following we shall argue that the theory 
is safe from tadpole up to dimension-eight level. 

First of all, it is not difficult to find out that under the symmetry Z^ 1 x Zf the allowed 
high- dimensional terms with S and Higgs bilinear term H u Hd take forms of S 2n+l (H U H d) m 
(with n and m being non-negative integers) which are odd powers of fields. So, according 
to the above discussions, they are safe operators. Next, let us consider the coupling of S 
with (N,N). There are two kinds of operators: one type are in the forms of S 2n+l (NN) m , 
which are safe; the other type are s 2n (N 2 ) 2m+1 (N 2 ) 2k (exchanging iV with N will not give 
anything new), which are potentially dangerous. But following the method of calculation 
given in [24], we can explicitly check that up to dimension-eight level (the next order is 
dimension-ten), the dangerous tadpole is absent. In fact, the inclusive potentially harmful 
terms are 



Then, due to the fact that in the super-Feynman diagram calculation one cannot connect 
two lines coming from the same vertex to form a loop (the superfield propagator which is 
proportion to 6(81 — 82) vanish), we conclude that they do not generate any tadpoles. For ex- 
ample, Fig. 1 will be the potential contribution to tadpole coming from the high- dimensional 
operator S 2 N 2 , but vanishes. Although we cannot ensure the safety of the theory up to ar- 
bitrary high orders, we believe that the sufficiently suppressed higher dimension operators 
will not destabilize the electroweak scale. 



S,L,Q^-(S,L,Q), 
H d ,N,U^-i(H d ,N, U), 
H u ,N,E c ,D^i(H u ,N,E c ,D). 



(18) 



S 2 N 2 , S 2 N 2 , S 2 N\ S 6 N 2 , S 2 N 6 , S 2 N 2 N 4 . 



(19) 




FIG. 1: The vanishing tadpole from operator S 2 N 2 . 
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III. PHENOMENOLOGICAL CONSEQUENCES 



In our model some new light degrees of freedom come into play around the weak scale and 
thus may cause some interesting phenomenology In this section, we give a preliminary study 
for such phenomenology We will mainly discuss those phenomenological aspects which are 
relevant to the recent experimental results. Some similar analysis has been carried out, 
albeit based on different models [26]. 



A. Light Sneutrino Dark Matter 



It is well known that in the MSSM the sneutrino has sizable coupling with Z— boson, 
which renders it not a good candidate for the cold DM [27] . To make it viable, several ways 
have been tried. One is to include lepton flavor violation and consequently the sneutrino 
is real and does not couple to Z diagonally [28]. Introducing right-handed neutrinos can 
also make the mixed v a good thermal DM candidate [29] (note that the analysis in [3] 
found that the almost pure right-handed sneutrinos can be thermal WIMP in the low scale 
seesaw, where the right-handed neutrinos have extra coupling to the NMSSM singlet S in 
a renormalziable way SN 2 and this can not only explain the TeV Majorana mass scale but 
also open effective annihilation channels for the right-handed sneutrinos). Another way is 
to abandon the canonical seesaw, and instead take the inverse seesaw to make v a viable 
DM candidate. A previous study (albeit not inspired by the recent detection experiments) 
has been carried out in [26]. In our following study we will show that in our model the v 
LSP can easily explain the DAMA/CoGeNT signals. Further, since the v LSP has distinct 
leptonic property, it may serve as a very special DM candidate favored by the recent indirect 
detection of DM. 

First, we check the slepton mass and pay special attention to a light DM candidate 
motivated by DAMA/CoGeNT. In the basis $ f = (v* L , N, N ), from Eq. (Al) we obtain the 
mass square matrix 



M~ « 



/ m\ + \M"z cos2/3 + m 2 D (— ttldAn + imtud cot (5) —tjidMn \ 

m 2 N + M 2 N + m 2 D B N 
\ rn% + M 2 N ) 



(20) 



where = AiM^ + nMff/Xi — AiA v 2 sin 2(3/2. For simplicity, we have assumed all pa- 
rameters to be real. The contribution from the small Majorana mass term for N and the 
corresponding soft bilinear term BfjN 2 can be safely ignored for almost all of the study 
on DM. However, it will induce a small but non-zero mass splitting for the CP-eigenstates, 
typically at the order B^/m^ ~ 10" 11 GeV, provided that a small bilinear parameter from 
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ordinary SUSY-breaking = XiAqv^/Mpi. This has important implication on the DM 
detection and will be discussed later. 




FIG. 2: The first two diagrams are the v\ — N scattering. The third diagram is the t— channel 
annihilation of v\ into a pair of r leptons. The last diagram is the t— channel annihilation of v\ 
into a pair of neutrinos via exchanging a neutralino (suppressed by neutralino mass). 

Now we focus on a light LSP ~ 8 GeV. We denote its fraction in various interaction 
eigenstates as 

v L = c L xV\ H , N* = c NX vx H , N = c^i H , (21) 

where dots denote other irrelevant heavy components. Typically, our desired LSP should 
have the doublet component c^i < 0.2 so as to explain the DAMA/CoGeNT signals by the 
Z— boson mediated large scattering rate o~^ p ~ 10~ 40 cm 2 [13] (note that the latest CDMS 
results seem to disfavor the DAMA/CoGeNT region [30]) and meanwhile give an effective 
coupling between the LSP and Z— boson gz^ZDivl with gz n ~ fi^cfj < 0.023 to satisfy the 
constraint from the Z— boson invisible decay width [31]. 

Two comments are in order. First, it is important that the LSP is complex so the diagonal 
coupling to the Z— boson is preserved. Consequently, the Z— boson just mediates a large 
elastic cross section uv l , P - This makes it distinguishable from other real or inelastic sneutrino 
LSP candidates. Second, although a large DM-nucleon scattering rate can be achieved via 
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Z— boson mediation, in principle the Higgs bosons such as h can also act as the mediator: 



where g qq h = h^xj cos/3 for down-type quarks and g qq h = h ul / sin/3 for up- type quarks. In 
the estimation, we assumed a large tan/3 = 30 and fc a = 0.14. Note that if the DM is 
light, compared to a fermionic DM (as heavy as v\) scattering with the nucleon, a scalar 
DM-nucleon scattering cross section is potentially enhanced by the coupling in the case of 
Z^/iii > A/fem^, where Xfh is the Yukawa coupling between the Higgs and fermionic DM. 
However, in this case the DM annihilates to the SM fermions like bb quite effectively (without 
p— wave suppressing, see Eq. (25)). As a result, it induces antiproton excess in the low energy 
spectra and may be disfavored by the PAMELA data [32]. 

We now discuss the way to give such a light sneutrino LSP with a doublet fraction 
cli ~ 0.2. If the SUSY-breaking is mediated by gravity, then all of the matrix elements 
are expected to be around the weak scale. Moreover, for the favored parameter pattern we 
have Mn ~ (9(500) GeV, so the scalar singlets N and N both have rather large mass term. 
As a result, the generation of a light sparticle needs some fine-tuning (reducing the mass 
eigenvalue by some proper non-diagonal elements). There are quite a few free parameters to 
tune so as to produce a light LSP required by DAMA/CoGeNT. Based on the structure of 
the mass matrix in Eq. (20), the following observations are useful to understand the origin 
of such a light LSP: 

• The complex LSP generically has significant z/£ component because various states are 
well mixed with each other. In practice, since is larger than and mo, in general 
vl is the dominant component of v\. Increasing m 2 L helps to increase the eigenvalue 
of the LSP and decrease the ul component. 

• To reduce the ui component, the element B N should be properly large to reduce the 
eigenvalues of the (2,3)-block so that v\ can mainly come from N and N. A large 
implies a large Ai or/and k/\\. But the element — mj^A^ + [imp cot (5 ~ —m^A^ is 
also important to determine the mixing. We assume A^ = Ai and k ~ A < 1. By 
the way, a small k may be required to generate a multi-TeV v s . 

• Although the soft masses m 2 N ,m 2 ^ Mjy, it is important to tune them to obtain a 
small mass of T>\ without significantly affecting cl\. 



.SI 

p,n 




( Phil \ 2 f(p,n) 2 ^ A'P> n ) n 

\m P , n + mvJ ^ fT « 9qqh+ 27 ^J T ° 9qqh 

^ \_q=u,d,s q=c,b,t 




(22) 
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Now we give a sample point: tan/3 = 3.0, A = 0.03, k = 0.08, i/n = 0.848, Ai = 0.14, 
A N = A l = 341 GeV, m N = m R = 171.2 GeV, m L = 400GeV, \i = 150 GeV. Other soft 
parameters such as squark masses are important to Higgs system but not important to the 
sneutrino sector, and thus they are fixed as 400 GeV. Then we get a snuetrino LSP: 

m Vl « 8.0 GeV, c L1 « 0.24. (23) 

Such a LSP does not lead to antiproton excess because its annihilation is mediated by 
Z— boson and is p— wave suppressed (see Eq. (24)). 

Now we check the annihilation rate of the sneutrino LSP pair, which determines the 
freeze-out relic density of LSP and its indirect signals. As a convention, we label v\ as 1 and 
its antiparticle as 2. Then we have to distinguish the annihilation modes: (anti)particle- 
(anti)particle self-annihilation rate o\\ = 022, and the antiparticle-particle annihilation rate 
C12 — a 2i- I n our interested parameter space, the contribution to cr 12 from s— channel 
Z— boson mediation is suppressed by (cli) 4 < 10~ 3 and also is p— wave suppressed, which is 
negligible. Explicitly we have 

7 771 ~ 

(* Z v rel ) ^(<W^M 2 (1 - 4a;,) 1 / 2 . (24) 

The contribution from the s— channel Higgs mediation is chirality suppressed. For example, 
the annihilation to bb is given by 

(WreiU - ^^h-A (1 " ml/ml) 3/2 
47r v z cos z p m\ 

~ 7.7 x ID- 6 x (^-Y ( cm 3 /*, (25) 

V 1.0 y V5GeV/ \ m h J 1 ' K ' 

where we fixed tan/3 = 30. If we take tan/3 = 3, then (a s v re i) b i will be reduced by about 
two orders. Thus, the antiproton production from DM annihilation is suppressed. 

The sneutrino LSP may have the leptonic annihilation property due to the significant N 
component and the survived vertex unLH u N at low energy. This vertex tends to play a 
major role in the annihilation of LSP for ~ 1. Explicitly, the average rate of annihilation 
into lepton pairs via t— channel chargino mediation is given by (the amplitude square is 
given in Eq. (A5)) 

(a t v rel ) fF ~ (^±^i v 2 + 4x/(x . _ Xf yj (1 _ 4x/) i/ 2; (26) 



where in the CM frame v re i = 2^1 — Ami~Js and 1 — vj = Axf ~ m^/m- ~ 10 1 for the 
GeV scale lepton. It is suppressed by an explicit chiral factor or p— wave suppression, but 
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still it dominates over the Higgs exchanging annihilation rate given in Eq. (25) in case of a 
light chargino Mq ~ rrih and/or a rather small fihii as well as a small tan/3. 

The annihilation rate presented in Eq. (27) has some implications. On the one hand, 
when we consider the DM annihilation today, DM is at the non-relativistic limit and thus 
it is safe to drop the terms oc v 2 ~ 1CT 6 . Then the chiral suppressed mode is barely able to 
account for the FGST: 

(27) 

On the other hand, at the DM decoupling era, the DM is semi-relativistic with v ~ 1/4 (for 
a typical WIMP m^/Tf ~ 20 — 25). Consequently, the p— wave annihilation mode is active 
and it is dominant for light leptons e and \i. So even the annihilation channel mediated by 
Higgs (to bb) is suppressed for the sake of suppressing antiprotons, we can still allow DM to 
have significant annihilation rate to the states such as ee (required by ADM discussed later), 
in additional to the fr channel with annihilation rate fixed by FGST. In other words, our 
model can realize the ambitious scenario that, the sneutrino LSP is an ADM thus predicted 
to be light, which supports both DAMA/CoGeNT and FGST. 

By the way, the rate an is from t— channel neutralino exchanging. It is suppressed by 
small mixing factor from (cjvi sin #i) 4 and thus the annihilation to light neutrinos is negligible 
in our model. 

B. Sneutrino as Asymmetric Dark Matter 

The sneutrino LSP can serve as the asymmetric dark matter (ADM). ADM is the most 
attractive framework to understand the coincidence between the density of the baryonic 
matter and dark matter, VL DM : Q b ~ 5 : 1 [16]. In some recent studies [17], usually one has 
to introduce extra DM candidate which interacts with leptons or quarks at high order so as 
to extend the lepton number or baryon number to the dark sector and moreover maintain 
the chemical equilibrium between the two sectors. DM and baryon number densities usually 
are roughly equal by virtue of the chemical equilibrium. Turn to our model, the sneutrino 
LSP carries lepton number and thus in principle it can serve as an ADM candidate (a 
previous study on such a LSP can be found in [18], which is based on the canonical seesaw 
equipped with special SUSY-breaking). This is a unique property of the sneutrino LSP 
because the other LSP candidates, the neutralino or gravitino DM does not carry lepton or 
baryon number. 

In our scenario we can assume the matter-antimatter asymmetry has been generated 
through some mechanism such as GUT-baryogenesis, Affleck-Dine mechanism or even EW- 
baryogenesis at the early universe. Especially, the lepton number asymmetry has been 
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generated and some is transferred to the (s)RHN system (N,N). The (s)RHNs decouple 
around Tf ~ <C T sp h, the out-of-equilibrium temperature of electroweak sphaleron 

process. From then on the baryon and lepton symmetry conserves respectively In such 
a way, the asymmetry stored in the sneutrino sector is calculable according to the well 
established chemical equilibrium condition. 

Obviously, it is crucial to determine the chemical equilibrium condition at T sp h- To do 
this, we follow the way developed in [33]. First, assign a chemical potential fj,^ to various 
particle <fi in our model. Then at temperature T, the asymmetry of particle <fi in thermal 
equilibrium can be expressed with (lower index is ignored) [34] 



where the Boltzmann suppressed factor fojim/T) denotes the threshold effect for heavy 
particle in the plasma. If the particle are in the ultra-relativistic limit m <C T, it tends to 
1; in the opposite limit, a decoupled particle has vanishing asymmetry density as expected. 
In this work, for simplicity, we take all sparticles, except for three lightest complex sleptons, 
to be highly Boltzmann suppressed and thus do not contribute to total charge asymmetry. 
In other words, as a quite crude approximation, at T sp h only SM and RHNs fermions as well 
as three light v particles store the charge asymmetry. 

Then the chemical equilibrium establishes relations between /x^. Around T sp h, at first all 
Yukawa couplings (including CKM mixing and dimensionless coupling in this model) have 
entered chemical equilibrium. As a result, the Higgs mediated flavor changing processes force 
all families to carry the same chemical potential. Next, the jy-boson and the charged Higgs 
mediated processes force the up and down components of the SU (2) l doublets to share the 
same potential (due to from those constraints, the quark/squark sector carries a potential 
fi UL while the lepton/slepton sector including all (s)RHNs carry fJL„ L ). Third, the gauginos 
mediate equilibrium for SM particle and its superpartner, while the soft trilinear terms make 
all sfermion families carry the same potential. Finally, the EW-sphaleron process is effective 
and thus the left-handed quarks and leptons satisfy 



Totally, there are two independent chemical potentials, one is fi VL and the other is the neutral 
Higgs /io- Especially, we have fj,„ L — — ^u L - Then the total baryon number is calculated 
to be 




f b (m/T) x fif (£) , (for bosons) 
f f (m/T) X 2|i (J ) , (for fermions) 



(28) 



0. 



(29) 



B ~ 



(30) 
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while the sneutrino asymmetry is given by 



2T 2 * 2kT 2 
i=l 

with z summing over all sneutrino states lighter than T sp (since all of their asymmetry will 
be finally transferred to LSP). Using Eq. (29) we get their net charge ratio 

n Zi _ S _ k . . 

which determines the relic number density ratio. Then from the observed energy ratio, we 
obtain 

m ~ = l*^L mp ~ 3.3 GeV - 10 GeV. (33) 

Taking the Boltzmann factor into account, one may get the right mass value of DM (~ 8 
GeV) required by the DAMA/CoGeNT and FGST. 

Of course, if the v is the viable ADM, its symmetric part should almost annihilate away 
during its freeze-out through v\ + v\ — > ff whose cross section is denoted as eri2, while 
the particle-particle annihilation rate an = <j<ii should be small (moreover, the particle- 
antiparticle oscillation by virtue of the dark lepton number violating term is negligible at 
this era, which will be addressed later). The only mode proceeds via exchanging a neutralino 
as stated before, which is highly suppressed. This is also consistent with a large Yukawa 
coupling Un. Note that this does not violate the FGST requirement since it only requires 
r ± -pair mode dominates over the qq modes while allows even a larger annihilation rate to 
the neutrinos as well as electrons (they do not produce hard gamma spectrum). 

One problem with V\ as the ADM may be that presently it may have no pair annihilation 
modes and thus has nothing to do with the gamma-ray signals. Amazingly, recall that v\ 
and v\ are just approximate degenerate states and they are split by the small sneutrino mass 
term coming from the Planck-scale-suppressed soft mass term. In the CP-basis, v\ and T>\ 
oscillate with a rate ~ 8m±. It just regenerates the symmetric DM after the LSP decouples. 
Note that the oscillating time scale t os ~ 5m^ should be much larger than the Hubble time 
scale tff ~ H~ l during the LSP freeze-out by symmetric annihilation (T12. Otherwise, the 
asymmetric part will be erased by significant oscillation effect, rendering the symmetric relic 
density determined by ordinary freeze-out dynamics again. In conclusion, the ADM sets a 
nontrivial constraint 

H(T = Tf) = I f < 5m ± (34) 
Mpi 

where 5m± = B^/m-p and z\ = m^/Tos ^> 1 are used. Roughly, it requires the soft term of 
S 2 N 2 /Mp\ suppressed by two orders compared to ordinary soft parameters. 
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C. Impact on Higgs physics 



In our model the Higgs phenomenology may be quite peculiar. First of all, v s is required 
to be several TeV in order to produce a proper neutrino mass scale (the NMSSM is near the 
decoupling limit). Given that \x = Xv s ~ 100 GeV, we immediately have A 1 and thus the 
lightest CP-even Higgs boson mass will not obtain enough enhancement from the tree-level 
quartic term via the coupling \SH u Hd- Actually, in the NMSSM the lightest CP-even Higgs 
boson mass square is approximately given by 

ml ~ Mf cos 2 2/3 + A V sin 2 2(3 + log f ) , (35) 

4ir z v z \ m t I 

where the last term is the dominant radiative correction (in the MSSM this term lifts up m\ 
to meet the LEP lower bound and thus needs a heavy stop quark mj ~ 1 TeV which causes 
the little fine-tuning problem). In principle, in our model a new contribution Am^ oc 
y 2 N m 2 D \og{m 2 ~ / m 2 D ) / IQit 2 from the coupling y^LH u N may be significant in some part of 
parameter space (e.g., y^ ~ 1) and thus help to alleviate the little fine-tuning problem 
(however, in the MSSM with the canonical seesaw, such a large Yukawa coupling of neutrinos 
just produces a negative effect and lowers the m h value [35]). 

But in this paper we focus on another more interesting effect on the Higgs, that may 
alleviate the little hierarchy problem. In the previous section, we have demonstrated the 
possibility of accommodating a rather light sneutrino LSP in the model, although at the 
price of some kind of fine tuning. In this section, starting from light sneutrino states (but 
not confined in the DM phenomenology), we pay special attention on the Higgs physics. 
Then, in the presence of SM singlets interacting with the Higgs sector, the decay modes of 
Higgs will be affected substantially if some states are sufficiently light. So the invisible Higgs 
decay scenario can relax the LEP2 lower bound on the Higgs mass [36, 37]. Notice that the 
sneutrinos are R— parity odd, their dominant decay mode is T>i — > LSP + X, where X is a 
neutrino for neutralino/gravitino LSP and a light Higgs (if kinetically allowed) or a pair of 
fermions for sneutrino LSP (this is natural when family structure is introduced). 

We denote the lightest Higgs component contained in the three CP-even Higgs as Rh u — 
h u \h + ■ ■ • (similar for Hd and S). Then the decay mode to bb has a width given by 



^ 66 ~^ d Vcos 2 /316^' 



N c = 3 is the color factor. This tree level value will be significantly reduced by the (super- 
symmetric) QCD correction at the NLO order. This decay mode is the dominant one in the 
MSSM. But in our model, new invisible decay modes are open: h — > ViVj and h — > v\i>\. 
The former is suppressed by small y 2 N and mixing factor, but the latter generically can be 
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rather significant. The decay to a pair of kinetically allowed sneutrinos has a width given 



where the coupling fi\ u is defined in Eq. (A2). From Eq. (A2) we see that this mode is 
potentially important, especially when the decays to neutrinos are suppressed by h 2 uX <C 1. 
Compared to the tree level bb decay mode in Eq. (36), this mode is enhanced by a factor 
{^hii/hdiyb rn 'h) 2 ■ Even without taking QCD correction into account, when tan/3 is in the 
small region tan /3 ~ 3 the Higgs invisible decay can be dominant decay mode. 

In light of the above general analysis, to show the possibility for the Higgs boson h to have 
a sizable invisible decay we consider the parameter space with the following characteristics: 

• To reduce the number of free parameters, all the new soft terms are assumed to be zero 
unless specified. This is well justified if the soft terms origin from gauge mediation 
(GMSB) or anomaly mediation (AMSB), where the soft terms involving singlets are 
proportional to Yukawa couplings of the singlets (in GMSB they are generated only 
through RGE effect and thus suppressed further by loop factor). Anyway, in the 
interested region 100 GeV ~mo« Mn, even if soft terms are generated at the weak 
scale, they do not play major role in our discussion. 

• At least one sneutrino T>i is lighter than m^/2 and has significant coupling //^/iz/jZ^. 
For simplicity, we use all the calculation in the previous section, and denote z/$ as v\. 
V\ is not necessarily the DM candidate and consequently many constraints are relaxed, 
such as the fraction of the i>l component in v\. 

• To keep /i^n significantly large, saying ~ 5 GeV, then from Eq. (A2) it can be clearly 
seen that mo and Mjv can not be too small. Generically, the lightest Higgs is h^R 
dominated, and the NMSSM decoupling limit shows a hierarchy h^i h u \ ^> h s ±. 
Then there are two obvious cases with physical simplicity. One is characterized by the 
^-dominated V\, but usually the singlet components in V\ are still appreciable and 
then the massive coupling is approximated to be 



which is rather large and sensitive to /i. This case can be realized clearly with a 
somewhat small ~ 150 GeV, properly reduced by the small mixing with heavy 
singlets. But > Mz/2 is required to satisfy Z— decay width. Oppositely, ~ 300 
GeV is larger compared to singlet mass term, leading to singlets dominated V\ and 



by 




(37) 



A*mi ~ c Ll m D (h ul m D /v + c m h dl ^/v) , 



(38) 
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H>hii ~ c^crnhdifnofi/v. By keeping cli < 0.25 and /i ~ 250 GeV, it takes value 
around 10 GeV. In this case the does not have to be tuned to set in the narrow 
window (Mz/2,rrih/2). But since mp is smaller than the previous case, it is expected 
that the Higgs invisible decay width will be suppressed more or less compared to small 
rriL scenario. 

Define R = T^^* /T h ^i b for a typical Higgs with mass m h ~ 100 GeV, where the 
decay rate is using the numerical result calculated in the NMSSMtools [38]. For the 
two cases, we consider the R varies with tan /3 with input parameters respectively 
given by: on the left panel of Fig. 3 we set mi = 160GeV, A = 0.024, k = 0.0593, 
and y N = 0.5, A x = 0.06; on the right panel, m L = 320GeV, A = 0.024, k = 0.0729, 
and Dm = 0.3, Ai = 0.0866. From the figure it is seen R decreases as (in larger tan/3 
region) tan (3 increases, and the Higgs invisible decay width benefits from a larger \i. 




FIG. 3: The lightest CP-even Higss boson h decays to a pair of light sneutrino with a ratio 
R defined as R = T^y^* /T h ^i b . From top to bottom, the curves corresponds to fixed /j, = 
220, 190, 160, 130, 100 GeV, respectively. 

• There is a third interesting case related with sneutrino DM. Typically, < Mjy, but 
we can still get a singlet-dominated T>\ by some proper tuning ( please see the detailed 
discussion in Section III A). In this dark matter inspired parameter scenario, both mo 
and take large values, leading to a rather large fihii- It implies that the Higgs 
inevitably decays to such LSPs with overwhelming branching ratio, typically > 90%. 
However, it should be stressed that to suppress the sneutrino LSP annihilation into 
bb, a small tan is favored. 

In conclusion, once there is light sneutrino states from this sector, we have to investigate 
their probably significant effect on the Higgs decay and their detection. By the way, we 
point out that a recent paper [39] studied an interesting phenomenology of a light CP-odd 
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Higgs and found it can decay to light neutrinos so as to relax the collider constraints on this 
particle. 

IV. CONCLUSION 

In this work, we dynamically realize the supersymmetric inverse seesaw mechanism in 
the NMSSM with a minimal field content. The crucial observation is that a x Z^ 1 sym- 
metry allows for a Planck-scale-suppressed lepton number violating dimension-five operator 
N 2 S 2 j Mpi, which can just produce a light neutrino mass scale m v ~ 0.01 eV. The theory 
accommodates the NMSSM ^-breaking domain-wall problem solution and is safe up to 
dimension-eight level. Naturalness and neutrino mass constrained the NMSSM parameters 
as v s ~ 5 x 10 3 GeV and A ~ 0.03. We then preliminarily explored some phenomenology of 
this model. In additional to the neutrino physics, the most salient thing is that the left-right 
mixed (with doublet slepton component en ~ 0.2) light sneutrinos with mass ~ 8 GeV can 
be the LSP candidate. Such a LSP dark matter can readily explain the DAMA/CoGeNT 
events and even explain the FGST gamma-ray signals. 

We also discussed the possibility that the sneutrino LSP is an asymmetric dark matter 
(ADM) and thus the ratio ft dm '■ ~ 5 : 1 naturally predicts a light DM. Further, the 
Higgs physics is also dramatically changed in our model, owing to the significant low energy 
coupling y N LH u N and the appearance of light sneutrinos. We found that the lightest 
CP-even Higgs boson probably dominantly decay to a pair of light sneutrinos and thus 
substantially relax the LEP bound. It provides a simple way to solve little hierarchy problem. 
Especially, in the region with light sneutrino LSP required by DAMA/CoGeNT, the Higgs 
boson invisibly decay with a dominant branching ratio (while the branching ratio to bb 
typically is less than 10%), which will affect significantly the hunting for the Higgs boson at 
LHC. 

Note that in the inverse seesaw mechanism the lepton flavor violating process like \i — > 
may be enhanced [40]. But in our model it depends on the flavor structure of (yN)ijH u LiNj 
when three families are incorporated. In principle, we can work in a limit {yN)ij ~ $ij to 
avoid such large lepton flavor violating processes. 
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Appendix A: Interactions and Annihilation of Sneutrinos 



In this appendix, we will present some necessary interactions involving the sneutrinos, 
including the mass terms, the coupling to Higgs as well as some amplitude square of the 
sneutrino annihilation. 

Collect all the relevant F— term and soft terms for the snuretinos, and write the doublets 



as Hi = (H+,H° U ),HJ = (H° d ,H~),L T 
the sneutrino system reads 



then the relevant scalar potential for 



V = \XxNN + kS 2 - \H° u H Q d \ 2 + \\SH% + y N v L N\ 2 
+ | - y N H° u u L + A^ivf + lAiSAf + \y N H° u N\ 2 

+ m 2 L \D L \ 2 + m%\N\ 2 + m%\ti\ 2 + (-y N A N H° u v L N + \ 1 A 1 SNft + h.c)j , (Al) 

where the D— term for L is not written. From the above we can read the mass matrix for 
the sneutrinos. We pay attention to the coupling between the lightest Higgs h and u^ui, 
namely the vertex ^///ai^i^i (l/\/2 factor stems from the kinetic normalization for the 
real Higgs h) with the coupling given by 



Phil =c N i/j,cos/3 



CLlh s 



m D 
v s sin f3 



VMN (h * R h > 

c ni— — {hdit&np - hui) 



+ Ki- 



rn 



D 



LI 



+ C 



A'l 



CLlCNl 



K 
A N 

m D m D h u i ' m D c N1 



fi h d i M N Cjv! 



+ h sl M_ 



N 



M 



N 



I 1 

[C m + CjyJ + C^Cni — 



cni v s 



2k 



(A2) 



Here the coefficients such as h u i are the /i-component contained in the real part of neutral 
Higgs R u . This coupling controls the strength of the LSP interacting with SM fermions via 
the Higgs h and plays a significant role in direct and indirect detection of DM. 

Next we calculate some amplitude square for V\ annihilation. First of all, we give the 
contribution from the t— channel chargino exchange, via a vertex in Eq. (8), 



Cy x D y N c Nl v?CP L T + h.c, 



(A3) 



with (C~ 



(H-,(H+y) and 



T 



(r L , (rjj)t) denoting the chargino and tau lepton, 



respectively. Here we neglected the contribution from the vertex g-iL *WL + c.c. because the 
light chargino is dominated by Higgsino and it also suppressed by smaller V\ from vl- In 
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the center mass frame, we denote the kinetics of initial DM pair and final fermion pairs as 



Pl =E(lJ), p 2 = E(l,-(3), 
k x = E?{\J'\ k 2 = E' (1,-/3'). 



(A4) 



For the particle % with mass m*, the velocity is defined as fa = \/l — 4m 2 / s with s 
(ki + k 2 ) 2 . Sum over the final spin states, the amplitude square is given by 

c2 



4M£ 



[8xf(xp 1 — xj) + v 2 (l — vj cos 2 9) — AvvjXf cos#] 



(A5) 



where the t— channel propagator factor is approximated to a contact form, and the terms 
linear in v will disappear after integrating over the angle to get the cross section. 

For the Z— mediated process, from the vertex c 2 L1 g 2 (yid^yl — v^d^y^Z^ , we have the fol- 
lowing p— wave suppressed (understood from the momentum dependent vertex proportional 
to (pi — P2) ~ v) amplitude square 

4 



M 



Ml 



2(pi -P2) ■ h(pi -p 2 ) ■ h - -(pi -p 2 f 



(Pi • kif + s\fi 



s 2 \v\ 2 



(1 + 2 cos 2 e) 



(A6) 
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